Bivariate Spline Method for Numerical Solution of Steady State Navier-Stokes Equations over Polygons in Stream Function Formulation

We use the bivariate spline nite elements to numerically solve the steady state NavierStokes equations. The bivariate spline nite element space we use in this paper is the space of splines of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the steady state Navier-Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth order equation, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in H2( ) of the nonlinear fourth order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C1 cubic splines is implemented in MATLAB. Our numerical experiments show that the method is e ective and e cient.